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Abstract. - The ongoing discussion whether thermodynamic properties can be extracted from 
a (possibly approximate) quantum mechanical time evolution using time averages is fed with 
an instructive example. It is shown for the harmonic oscillator how the Hilbert space or an 
appropriately defined phase space must be populated in terms of coherent states in order to 
obtain the quantum result respectively the classical one. 



Introduction. — 

Transport models which employ wave packet dynamics are widely used in atomic, solid state 
and nuclear physics today. Examples range from the description of hydrogen plasmas |^ over 
the investigation of critical properties of noble gases ||] to simulations of heavy ion collisions 
in nuclear physics ^, |^. In all of these models the quantum mechanical time evolution 
is approximated by the evolution of a (possibly antisymmetrized) product state composed of 
Gaussian single-particle wave packets (squeezed coherent states ^) 

/ ^1 ^ ^ \ f * ^ ^1 

(x\r,p,a} cx exp<^ — + -p ■ x > , (1) 

where the time dependence of the state is given through the time dependence of the parameters, 
as there are the mean position r , the mean momentum p and in some cases the complex width 
a. 

Although these models were designed to describe non-equilibrium situations their thermody- 
namic equilibrium properties are important for two reasons. Firstly, their statistical properties 
have great influence on final observables like fragment distributions in collisions. Secondly, 
long time evolutions are used to determine equilibrium behavior and to estimate quantities 
like the degree of ionization of a plasma [Q or the caloric curve for the nuclear liquid-gas phase 
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transition 0|. The last point is of special interest since approximate quantum evolutions are 
a possible way to examine thermodynamics properties in cases where the energy spectrum 
and thus the partition function is not known. In these cases the time-dependent Schrodinger 
equation cannot be solved, too. 

Whereas thermostated time evolutions are on firm grounds in classical mechanics ^, 10 



and have been used successfully as well for equilibrium, e.g. to investigate classical spin 
systems, as for non-equilibrium, e.g. to study glass transitions, in quantum mechanics less 
attempts were made, some of them are given in refs. 0, ^ ^ |3[ |^ |l5[ |l6[ ^ 
Following these articles one realizes that the matter is still under debate. The questions 
are twofold. The first concerns the problem whether it is in principle possible to extract 
thermodynamic properties from a quantum mechanical time evolution. The second question 
asks in how far a coherent state used as a trial state in an approximate dynamics like a quantum 
molecular dynamics is able to visit the Hilbert space in time with the correct weight. 

Some of the areuments, which doubt a success of time averaging are loosely speaking like 
(see for instance ^ , 20 ) : Since the equations of motion, either approximative ones like 



in Antisymmetrized Molecular Dynamics (AMD) or Fermionic Molecular Dynamics 

(FMD) [^/ or the Schrodinger equation itself, have a symplectic structure, like Hamilton's 
equation of motion in classical mechanics, a time averaging should lead to classical statistics at 
its best. The dynamics does not know about the discrete structure of the energy levels, it is not 
quantized. One can excite a system, which is described by coherent states, by any small energy, 
so it does not know about the finite spacing of energy levels. Coherent states are quasi classical 
states, the dynamics must lead to classical statistics. Or last but not least, the restriction to 
product states is equivalent to a mean field calculation, the result of time averaging is therefore 
classical. This ongoing discussion got new momentum during the program INT-98-2 [p2[ at 
the University of Washington in summer 1998. 

It is the aim of this letter to demonstrate with the algebraicly simple, but instructive 
example of the harmonic oscillator how the Hilbert space or an appropriately defined phase 
space must be populated in terms of coherent states in order to obtain the quantum result or 
the classical one. That means, that one cannot a priori say, that a (possibly approximate) 
quantum mechanical time evolution is not able to reproduce the correct ensemble average nor 
that the restriction of the trial state to wave packets or products of them impedes correct time 
averages. 

The harmonic oscillator is not a mere example, nowadays it is of special interest for the 
dynamics of bosonic or fermionic atoms contained in magnetic traps |2|, [2^, |2^ as well as 
for all systems that show an equidistant level spacing close to the ground state like nuclei or 
Luttinger liquids. 



Thermodynamic mean. - 

Given the Hamilton operator H of the harmonic oscillator 

H = huj I a^a + i ) , (2) 



2, 

coherent states are defined as eigenstates of the destruction operator a 



Each coherent state is characterized by a complex parameter z which corresponds to a pair 
of real parameters {r,p). Coherent states span the Hilbert space, they are over-complete and 
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obey the completeness relation 



1(1) = / \z){z\ , d2z = dRe(z) dlm(z) (4) 



TT 



A thermodynamic mean of an observable B in the canonical ensemble is given by the trace of 
this operator together with the statistical operator. This trace can be expressed in any basis, 
it needs not to be the eigenbasis of , it may be a basis like those of coherent states which 
are characterized by a continuous parameter |21, ^ 



1 f d 



Z{p)J TT 



— e-^ ^-\z?{i-e-f'--) ^ g-i I B I P^^z ) 



1 r A^z 



Z{(3) 

The partition function then reads 



m = / —{z\e-''i\z) (6) 

TT 
TT 

Given the excitation energy 

n{z) = {z\H\z)-^huj = huj\z\^ = ^ + ^mu;'r\ (7) 

the weight with which a coherent state \ z) = \r,p) contributes to a thermodynamic mean 
is just 

wquim = e-l^l^(^''^^-i)=e-«WK'^-i)/(^-) (8) 

where all terms depending only on the zero point energy have been omitted since they cancel 
with the respective terms in the partition function. 

One can now interpret the space of parameters z = (r, p) as a phase space and reformulate 
eqs. (H) and (^) 

ZW) = / i»qm(3) . (10) 

Then Wqn\{P) is the thermal weight in this phase space and B{r,p) = ( t'jP | B \ r,p) a function 
of the phase space variables. 
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Fig. 1. - Quantum weight (solid line, eq. (^)), classical weight (dashed hne, eq. (|ll[)), next order 
(dashed dotted hne, eq. (|l2[)) for the harmonic oscillator potential. 



Classical limit. - 

The connection to classical mechanics is established either by performing the classical 
{h ~* 0) or the high temperature limit (/? — > 0). The exponential in the exponent of Wqm(/^) 
can be expanded, huj drops out and the weight (|^) approaches the classical result 

w^liP) = e-^«(^) . (11) 

Keeping the next order in the Taylor expansion of the exponential function one gains 

WnoiP) = e-^(i+i ^^■^)^^'-^ . (12) 

Figure]^ shows all three weights for a specific inverse temperature /3 — l/{huj). One sees that at 
the same temperature the quantum mechanical distribution is much narrower in jzp oc Ti.{z) 
than the classical one, meaning that the quantum mean excitation energy is less than the 
classical mean energy. 

Conclusions. - 

The purpose of the above exercise is to remind that one can represent quantum mechanics 
and statistics in a basis which is characterized by a continuous parameter. Although this 
parameter does not know anything about the discrete nature of the energy eigenvalues and 
moreover leads to the temptation of a classical interpretation, the quantum statistical results 
are correct, because they do not depend on the chosen representation. 

Therefore, if the time evolution of a particle enclosed in a harmonic oscillator and described 
by a coherent state populates the Hilbert space, or phase space in the sense of (^), with the 
weight wqmiP) the statistical properties are quantum mechanically correct. 

It is a different question whether a specific approximation of the time-dependent Schrodinger 
equation or a specific coupling to a thermostat actually results in such a time evolution. In 
the spirit of Nose Jill one could for example define a thermostat which works by scaling 
coordinates and momenta so that the coherent state visits the Hilbert space according to the 
desired weight. One possible Nose Hamiltonian would be 



J.SCHNACK THERMODYNAMICS WITH COHERENT STATES 



5 



L 

3 


1 


1 


1 1 1 

(r^.Px.ry.Py) ... ■■ - 










1— 








3r-2 
LU 













1 


■ lr/P,a> 

1 1 1 " 




0.0 


0.5 


1.0 1.5 2.0 



Fig. 2. - Excitation energy for different degrees of freedom. The solid line shows the classical result 
for a one-dimensional harmonic oscillator, the dotted line the result for two-dimensions; the dashed 
line displays what happens if the coherent state is replaced by a wave packet with time-dependent 
width. 



Another example, which is formulated with the help of a thermostat for a constrained dynam- 
ical system, is given in ref. p2| . 

Concerning approximations of the time-dependent Schrodinger equation it could be shown 
numerically that distinguishable particles contained in a common harmonic oscillator field and 
interacting via a short-ranged two-body force which are described by coherent states approach 
classical statistics on time-average but that the extension of the coherent state towards 
wave packets with time-dependent width (squeezed coherent states) as used in FMD [El] 



(a;|r,p,a) cx exp | ^ j ^'^^^ 

is sufRcient to obtain the quantum result. Having a closer look on the excitation energy 

'H{r,p,a) = {r,p,a\ H \r,p,a) - ^huj (15) 
■m,uj r H — — I- muj — - — nuj 



2m 2 2 \2mRe{a) 2 Re (a) / 2 

which now depends not only on the degrees of freedom r and p, but also on the width 
parameter a, this might seem counter intuitive compared to classical statistics. As illustrated 
in fig. ^ in classical statistics additional degrees of freedom, for instance the extension of 
the one-dimensional oscillator to two dimensions, lead to higher mean energy at the same 
temperature, whereas the enlarged freedom of the wave function introduced by the non-classical 
width parameter does the opposite. The wave function is now able to "feel" the underlying 
density of states which is less than the classical one, therefore, less excitation energy belongs 
to the same temperature. 

In general one must say that whether an approximate quantum dynamics is able to reflect 
thermal averages does depend on the freedom of the trial state and also on the observables 
one is interested in. It may very well be that the parameterization of the trial state is rich 
enough for single-particle observables but not for many-particle observables. 
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